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Q ! Use of the highly sensitive Hokupa'a curvature wavefront sensor has allowed 

^ ! for the first time direct adaptive optics (AO) guiding on brown dwarfs and VLM 

^! stars {SpT = M7 — L2). An initial survey of 9 such objects discovered one 

O ; 0.15" binary (2MASSJ 1426316+155701). The companion is about half as bright 

^ ! as the primary {AK = 0.61 ± 0.05,Aif = 0.70 ± 0.05) and has even redder 

! colors H — K = 0.59 ± 0.14 than the primary. The blended spectrum of the 

! binary has been previously determined to be M9.0. We modeled a blend of an 

^ i M8.5 template and a L1-L3 template reproducing a M9.0 spectrum in the case 

. .c^ ; of AK = 0.61 ± 0.05,Aif = 0.70 ± 0.05. These spectral types also match the 

observed H-K colors of each star. Based the previously observed low space mo- 
tion and Ha activity we assign an age of 0.8~^QlGyr. Utilizing this age range and 
the latest DUSTY models of the Lyon group we assign a photometric distance 
of 18.8;j;^^pc and masses of Ma = O.O74l°:[;?^M0 and Mb = O.OQQto fiMQ. We 
therefore estimate a system separation of 2.92^q^IAU and a period of IS.S^l'Hyr 
respectively. Hence, 2M1426 is among the smallest separation brown dwarf bi- 
naries resolved to date. 



Subject headings: instrumentation: adaptive optics — binaries: general — stars: 
evolution — stars: formation — stars: individual (2MASSJ 1426316+155701) — 
stars: low-mass, brown dwarfs 
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1. Introduction 

Since the discovery of Gl 229B by Nakajima et al. (1995) there has been intense interest 
in the direct detection of brown dwarfs. To fully understand this new population of sub-stellar 
objects it is critical to characterize binary brown dwarfs (systems in which both components 
are likely substellar) . Brown dwarf binary systems are of critical importance since they alone 
allow the masses and luminosities of brown dwarfs to be measured directly (Kenworthy et 
al. (2001), Lane et al. (2001)). Without direct dynamical mass measurements theoretical 
evolutionary models of brown dwarfs (Chabrier et al. (2 001), Burrows et al. (2000)) cannot be 
calibrated. Therefore, there is great interest in directly detecting close brown dwarf binaries. 

Moreover, the binary frequency of brown dwarfs is interesting in its own right, since little 
is known about how common binary brown dwarf systems are. It is not clear currently if the 

brown dwarf binary distribution in mass ratio (q) and separation is similar to that of M stars; 
in fact, there is emerging evidence that double brown dwarfs (hereafter called binary brown 
dwarfs) tend to have mass ratios close to unity and small separations (Martin, Brandner, & 
Basri (1999)). However surveys also suggest the binary frequency itself is similar to that of 
M stars (Reid et al. (2001a)). 

Despite the strong interest in such binary brown dwarfs, very few such systems have 
been detected to date. Before HST/NICMOS ran out of cryogen the first double L dwarf was 
detected (Martin, Brandner, & Basri (1999)). A young spectroscopic binary brown dwarf 
(PPL 15) was detected in the Pleiades (Basri & Martin (1999)), but this spectroscopic system 
is too close to get separate luminosities for each companion. However, a large HST/NICMOS 
imaging survey by Martm et al. (2000) of VLM dwarfs in the Pleiades failed to detect any 
brown dwarf binaries with separations > 0.2" (>~ 27 AU). Detections of closer binary 
systems were more successful. The nearby object Gl 569B was resolved into a 0.1" (1 AU) 
binary brown dwarf at KECK and the 6.5m MMT (Martin et al. (1999); Kenworthy et 
al. (2001); Lane et al. (2001)). Keck seeing-limited NIR imaging marginally resolved two 
more binary L stars (Koerner et al. (1999)). A survey with WFPC2 detected 4 more (3 
newly discovered) close equal magnitude binaries out of a sample of 20 L dwarfs (Reid et al. 
(2001a)). Then Potter et al. (2001) guided on HD130948 with adaptive optics and discovered 
a binary brown dwarf companion. Hence, the total number of binary brown dwarfs known 
is just nine. In addition, of these, only eight have luminosities known for each component, 
since one is a spectroscopic binary. 

In this paper we present a newly discovered binary brown dwarf 2MASSJ 1426316+155701 
(hereafter 2M1426) which is a significant addition to this list since it is the only example 
of a M8.5 primary with a L1-L3 secondary which is a large mass ratio compared to other 
brown dwarf binaries. At the observed separation of 0.15", the 2M1426 binary is the 4*'* 
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smallest separation resolved brown dwarf binary currently known, and so it should likely 
play a significant role in the mass luminosity calibration for brown dwarfs. Moreover, we 
show for the first time that it possible to guide directly on a low mass star/brown dwarf 
with adaptive optics on an 8m telescope, opening up a new technique for diffraction-limited 
brown dwarf observations. 



2. Adaptive Optics Imaging with Gemini &; Hokupa'a 

Our approach to the detection of brown dwarf binaries is somewhat different from that 
tried before. We wished to utilize both diffraction-limited resolution and NIR (1 — 2.5/im) 
imaging at the wavelengths where brown dwarf spectral energy distributions peak. To detect 
the closest binaries we utilized the largest aperture possible (D~ 8 m). Hence ground-based 
NIR imaging with adaptive optics (AO) was a logical choice (cf. Close 2000 for a review 
on AO). However, due to the extreme optical faintness of these brown dwarfs {V ~ 20) 
most AO systems can not lock on such faint targets. The exception are curvature based AO 
systems which employ red-sensitive photon counting avalanche photodiodes (APDs) in their 
wavefront sensors. Such a sensor can lock on a target as faint as V=20 as long as it is very 
red (V-I=4). Since nearby VLM stars/brown dwarfs of spectral type M8-L1 have V ^ 20 
and y — / ~ 4 we decided such targets would be possible with a curvature AO system on 
an 8 meter telescope. 

Recently, the University of Hawaii curvatTirc AO system Hokupa'a (Graves et al 1998; 
Close et al. 1998) was moved to the Gemini north telescope and provided as a visitor 
instrument to the Gemini community. Although Hokupa'a's 36 elements were too few to get 
very high Strehl images on an 8m telescope, it was well suited to locking onto nearby, faint, 
red M8-L1 stars and producing 0.1" images (which is close to the 0.05" diffraction-limit at 
H band) . We decided to utihze this unique capability to survey the nearest extreme M and 
L stars to characterize the nearby binary brown dwarf population. 

Here we report the results of our first observing run on June 20, 2001. We targeted 
VLM stars and brown dwarfs identified by Gizis et al. (2000). In total, 9 brown dwarfs 
were observed over the first half of the night. Typically we achieved 0.1 — 0.15" images of 
~ 10 min total integration time in the H (1.6/im) band. Since the survey is on-going we will 
report in more detail about all the objects observed in a more comprehensive paper in the 
near future. In this paper we limit the discussion to the 2M1426 binary observations. 

One out of 9 of our targets (2M1426) was clearly a tight binary with separation ~ 0.15". 
We observed the object dithering over 4 different positions on the QUIRC NIR 1024x1024 
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detector with 0.0199"/pix (Hodapp et al. (1996)). At each position we took 3x5s exposures 
and 2x60s exposures. In addition, we took a series of lOxlOs K and K' images all in the 
same position. 

3. Reductions 

We have developed an AO data reduction pipehne in the IRAF language which maxi- 
mizes sensitivity and image resolution. The pipeline cross correlates each image, then rotates 
each image so north is up and east is left, then median combines the data with an average 
sigma clip rejection at the ±2.5a level. By use of a cubic-spline interpolator the script pre- 
serves image resolution to the < 0.02 pixel level. Next the script produces 2 final output 
images, one that combines all the images taken, and another where only the sharpest 50% 
of the images are combined. The final images have a FOV of 30x30". 

This pipeline produced a final deep 480s H image (8x60s; FWHM = 0.197"), a final 
unsaturated 60s H image (12x5s; FWHM = 0.131"; see Figure 1), and a unsaturated 50s K 
{FWHM = 0.197") and K' (lOxlOs; FWHM = 0.197") image. 

4. Analysis 

In Table 1 we present the analysis of the images taken of 2M1426. The photometry 
was based on DAOPHOT (Stetson (1987)) PSF fitting photometry. The PSF used was the 
12x5s unsaturated data from the next brown dwarf observed after 2M1426. This PSF "star" 
had a similar brightness (K=12.5) and late M8.5 spectral type and was observed at a similar 
airmass. The flux ratio measured by DAOPHOT was normalized by the total flux of the 
binary in a 15" aperture. The resulting magnitudes are listed in Table 1. 

The platescale and orientation of QUIRC was determined from a short exposure of the 
Trapezium cluster in Orion and compared to published positions as in Simon Close & Beck 
(1999). From these observations a platescale of 0.0199"/pix and an orientation of the Y-axis 
(0.3 degrees E of north) was determined. 

The astromctry was based on lightly (100 iterations) LUCY deconvolved data from 
the 12x5s unsaturated H image (see Figure 1). For the deconvolution we used the same 
PSF as above. A good deconvolution result was obtained with a restored resolution of 
FWHM — 0.080". This high resolution allowed for excellent astrometric accuracy to be 
achieved by centroiding on the well separated stars. These astrometric measurements were 
also checked against the PSF fitting photometry of DAOPHOT and were found to be in 
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excellent agreement (DAOPHOT gave separation = 0.155", PA = 343.5°). 

5. Discussion 

5.1. Is the companion physically related to the 2M1426 system? 

We believe there is a very high probability that the companion is physically associated 
with the 2M1426 system. This is likely due to the small space density of very red {H — K = 
0.59) background objects in the field. Furthermore, in the 9000 square arcsecs already 
surveyed we have not detected a similarly reddened background object in any of the fields. 
Therefore, we estimate the probability of a chance projection of such a red object within 
0.15" of the primary is < 1x10"^. We conclude that this very red, cool object is physically 
related to the 2M1426 primary and hereafter refer to it as 2M1426B. 

5.2. What Eire the spectral types of the components? 

We do not have spatially resolved spectra of both components; consequently we can only 
try to "blend" spectra of a hotter and cooler star while preserving the observed AK = 0.61 
and AH — 0.70. In Kenworthy et al. (2001) it was found that a M8.5 and a L1-L3 spectrum 
could be blended together (with AK — 0.61) to produce a M9 "blended" spectrum. Since 
Gizis et al. (2000) observed 2M1426 to have a blended M9.0 spectrum we can reasonably 
assume that the components may have spectral types of M8.5 and L1-L3. Hence we estimate 
(until a proper spatially resolved spectrum is obtained of both components) that 2M1426A 
is similar to a M8.5 and 2M1426B is similar to a L1-L3 spectrum. Furthermore, the H — K = 
0.47 color of the primary is consistent with a MS. 5 SpT. Similarly, the H — K = 0.59 color 
of the companion is consistent with a L1-L3 SpT Reid et al. (2001b). So we will assume a 
MS. 5 SpT for the primary and a L1-L3 SpT for the companion until proper spatially resolved 
spectra are obtained and the true spectral types determined. 

5.3. What is the distance to the 2M1426 system? 

Unfortunately there is no published parallax to the 2M1426 system. We can estimate, 
however, the distance based on the spectral types and a range of possible ages for the 
system. To do this, calibrated theoretical evolutionary tracks are required for objects in 
the temperature range 2300-1800 K. Recently such a calibration has been performed by 
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two groups using dynamical measurements of the G1569B brown dwarf binary. Prom the 
dynamical mass measurements of the GL569B binary brown dwarf (Kenworthy et al. (2001), 
and Lane et al. (2001)) it was found that the Chabrier et al. (2001) and Burrows et al. 
(2000) evolutionary models were in reasonably good agreement with observation. Based on 
the latest Chabrier ct al. (2001) DUSTY models (see Figure 2) we find that even if the age 
of 2M1426 was as old as 7.5Gyr the closest it could be is 17.8 pc based on the observed 
Ka = 12.16 and 5pT^=M8.5. On the other hand, if 2M14216 is as young as 0.5 Gyr it could 
be as far as 20.3 pc with Ka — 12.16 and S'pT^=M8.5. Since it is highly hkely that the age 
of 2M1426 is between 0.5-7.5 Gyr, we adopt a distance range of 17.8-20.3 pc. 



5.4. What is the age of the 2M1426 system? 

Estimating the age for the 2M1426 system is difficult since there are no Li measurements 
yet made. The age could be anywhere from 0.5-7.5 Gyr. However, the low proper motion 
observed by Gizis et al. (2000) for the system {/iRA = 0.108"/yr, i^DEC = -0.056'7yr) 
gives Vtan — 12 km/s. Unfortunately the low resolution of the spectrum obtained by Gizis 
et al. (2000) does not allow the radial velocity (and hence full 3-dimensional space motion) 
to be calculated at this time. However, such a low Vtan = 12 km/s velocity makes 2M1426 
the lowest velocity M8-M9 in the entire survey of Gizis et al. (2000). This strongly suggests 
a young age since the system has not developed a significant random velocity like the other 
older M8-M9 stars in the survey. Moreover, the low 1.2 Angstrom EW of Ha may also 
suggest a youthful age for such a late spectral type (where Ha appears more strongly in 
older systems; Gizis et al. (2000)). Indeed, 2M1426 occupies a uniquely young position as 
the lowest Ha and lowest Vtan of any of the 14 M8-M9 objects in the Gizis et al. (2000) 
survey. Therefore, we assign a young age of 0.81q;3 Gyr to this system, but leave large error 
bars. 



5.5. How does 2M1426 compare to the other binary dwarfs known? 

As was mentioned in the introduction there are 8 binary brown dwarfs (published at 
this time) with resolved components (i.e. non-spectroscopic) . A listing of these in order of 
separation can be seen in Table 3. 

Prom Table 3 we see that, other than the very close (lOpc) system GL569B, the next 4 
widest systems all have estimated periods in the range 11-16 yr. These periods in Table 3 are 
crudely estimated by assuming face-on circular orbits. However, even such crude estimates 
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do allow us to see that 2M1426 is in the top 5 of the shortest period brown dwarf binaries 
currently known. This is of great importance since it is only those systems with separations 
less than 4 AU that will have periods less than 25 yrs. Since a good fraction of a period is 
required to calculate the dynamical mass, such objects are of great importance. The shortest 
period resolved brown dwarf binary (Gl 569B) has already had its orbit calculated and the 
first zero-point to the theoretical tracks has been applied (Kenworthy et al. (2001); Lane et al. 
(2001)). However, one measurement does not cahbrate the 3-dimensional Mass-Luminosity- 
Age tracks (see Figure 2) , and so dynamical masses and luminosities are required for more of 
the systems in Table 3. Hence, the short period binary 2M1426 makes a significant addition 
to the known brown dwarf binaries. 



5.6. The masses of the components 

Unfortunately no Li observations exist yet for the system, so we must use the models in 
Figure 2 to estimate the masses. We can estimate the masses of the components based on the 
age range of 0.5-7.5 Gyr. The models suggest a mass range of Ma — 0.063 — 0.079Mo and 
Mb = 0.051 - O.O72M0 for our spectral types of M8.5 and LI and the age range of 0.5-7.5 
Gyr, respectively. At the older ages (> IGyr) the primary may be on the stellar /substellar 
boundary, but the substellar nature {Mb < O.O75M0) of the companion is certain even in 
the unlikely age 7.5 Gyr, according to the models. 

5.7. Future observations 

Future observations of this binary should determine if there is still Li present in its 
spectrum. The best Li observation would spatially resolve both components in the visible 
spectrum. Currently there is no visible wavelength AO systems capable of guiding on a 
~ 20 source and obtaining resolutions better than 0.1" (Close 2000). Therefore, we 
will have to carry out such observations from space with HST/STIS perhaps. Parallax 
measurements should also be obtained in the near future. Within the next few years one 
should also be able to measure the masses of both components as they complete a significant 
fraction of their orbit. Hence, further observations of 2M1426 are a high priority for the 
calibration of the brown dwarf mass-age-luminosity relation. 

We thank Isabelle Baraffe for supplying us with theoretical models computed by herself, 
France AUard, Gilles Chabrier and Peter Hauschildt. Mathew Kenworthy helped produce 
Figure 2. LMC acknowledges support by the AFOSR under F49620-00- 1-0294. Observa- 
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observations were supported by the University of Hawaii AO group (D. Potter, O. Guyon, 
P. Badouz and A. Stockton). Support for Hokupa'a comes from the National Science Foun- 
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Table 1. 


The 2M1426 System June 20, 2001 


Olwcrvcd 


Value 


Calculated 


Value 


PaiauKM 




Paraiii('l(n' 




AK 
AK' 
AH 

Ang. Sep. 
PA 


0.61 ± 0.05 
0.65 ±0.10 
0.70 ± 0.05 

0.152 ± 0.005" 
344.05 ± 0.5° 


Age 

Distance* 
Separation 

Mtot 

Est. Period 


0.8t«-^Gyr 
18.8ll;^^pc 

2mil-llAu 

0.1401;J;«i1Mq 



''photometric distances from the models of Chabrier et al. (2001) 
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Table 2. Summary of the 2M1426 A & B components 



Parameter 2M1426A 2M1426B 



H mag 


12.63 ± 0.05 


13.34 ± 0.10 


Ks mag 


12.20 ± 0.07 


12.80 ±0.14 


K mag 


12.16 ±0.05 


12.75 ±0.10 


H-K color 


0.47 ± 0.08 


0.59 ±0.14 


Est. SpT 


M8..5 


Ll-L.'? 


Mass 




o.o(i(ii;;;;;^^;Afe 



-13- 



Table 3. Published Resolved^Brown Dwarf Binaries 



Name 


Separation 


Est. Period'' 


Ref.<: 




All 






Gl 569Bd 


1.0 


3 


1,2 


HD130948B'i 


2.5 


11 


3 


2M0746 


2.7 


12 


4 


2M1426 


3.1-2.7"^ 


16-12'= 


this paper 


2M0920 


3.2 


16 


4 


2M0850 


4.4 


28 


4 


DENIS 1228 


4.9 


33 


5 


2M1146 


7.6 


63 


6 


DENIS 0205 


9.2 


72 


6 



''resolved binaries have measured fluxes for each compo- 
nent, the one spectroscopic binary PPL 15 (Basri &; Martin 
(1999)) is excluded 

'"This "period" is simply an estimate assuming a face-on 
circular orbit 

'=REFERENCES-(1) Kenworthy et al. (2001); (2) Lane 
et al. (2001); (3) Potter et al. (2001); (4)Reid et al. (2001a); 
(5) Martin, Brandner, & Basri (1999); (6) Koerner et al. 
(1999) 

'^Gl 569B and HD130948B are binary brown dwarfs that 
orbit normal Stars, for AO observations these bright pri- 
mary stars were guided on, not the brown dwarfs 

''values for age range 0.5-7.5 Gyr respectively 
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Fig. 1. — In figure (a) wc sec the 5s exposure 2MASSJ 1426316+155701 binary at an H band 
resolution of 0.131". In Figure (b) we show the system after being LUCY restored to 0.080" 
resolution. Only astrometry was derived from the deconvolved images, photometry was 
derived from DAOPHOT PSF fitting of the undeconvolved data. The pixels are 0.0199"/pix. 
The contours are linear at the 95, 75, 55, 35, 15% levels. 
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Fig. 2. — The latest Chabrier et al. (2001) DUSTY evolutionary models, and the temperature 
scale of Leggett et al. (2001). The locations of the 2 components of 2M1426 are indicated 
by the sohd circles. Note the good agreement of the secondary's position with the 0.8 Gyr 
isochrone and the 1900K (~L1) isotherm. The large error bars are due to the errors in the 
uncertain age of the system O.SlgigGyr. 



